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Abstract 

^ , We search for approximate, but analytic solutions of the pairing 

^ ' problem for one pair of nucleons in many levels of a potential well. 

C ■ For the collective energy a general formula, independent of the details 

of the single particle spectrum, is given in both the strong and weak 
coupling regimes. Next the displacements of the solutions trapped 



X 



^ • in between the single particle levels with respect to the unperturbed 



energies are explored: their dependence upon a suitably defined quan- 
tum number is found to undergo a transition between two different 
regimes. 

PACS: 24.10. Cn; 21.60.-n 
Keywords: Pairing interaction 

1 Introduction 

In this paper we deal with the problem of the pairing Hamiltonian for one 
pair of nucleons living in a set of levels of a potential well. 
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There are several motivations for this study. First, to find approximate, 
but analytic, solutions for the eigenvalues and the eigenfunctions, which we 
believe to be not only interesting per se, but also useful for paving the way 
to the general problem of n interacting pairs (see ref. [||). Second, to con- 
nect the energy of the collective mode, both in the strong and in the weak 
coupling regime, to global features of the single particle levels spectrum like 
the variance, the skewness, the kurtosis, etc. Third, to unravel the remark- 
able pattern displayed by the solutions trapped in between the single particle 
levels, already hinted at in ref. p. Indeed by connecting the trapped solu- 
tions to a quantum number A, it is found that their behaviour versus A is 
not only smooth, but displays a transition between two different regimes. 
Interestingly, this transition may be on one side related to a sort of sum rule 
obeyed by the trapped eigenvalues (stemming from the Viete conditions for 
the solutions of an algebraic equation) and on the other to the basic nature 
of the pairing interaction. 

The paper is organised as follows: in section ^ the general formalism 
is presented and the basic equations are deduced in the framework of the 
Grassmann variables; in section ^ the collective solution is addressed in both 
the strong and weak coupling regimes; in section ^ the trapped solutions 
and their regularities are studied within the harmonic oscillator and related 
models for the single particle potential well. 

2 General formalism 

Consider a shell of an average potential well, whose shape needs not to 
be specified, with L non-degenerate single particle levels of energy and 
angular momentum j^, the associated multiplicity being 2Qi, = + 1 
(1 < z/ < L). 

Let identical fermions (e.g., neutrons) living in this set of levels interact 
through the pairing force. The Hamiltonian of the system H = Hq + Hp 
splits then into a single particle 

HQ = Y^e^ 5Z 4.m,>>m. (1) 
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and into a pairing interaction 

L 

Hp = -Gj2 (2) 

fJ,,U=l 

term. In (0) 

Ojm, being the nucleon's destruction and creation operators. The oper- 
ators and AJ^ destroy and create pairs having total angular momentum 
J = in the level j^. 

As well-known, the problem of the Hamiltonian (|l],0) has been addressed 
by first diagonalising the associated bosonic Hamiltonian and then by 
accounting for the Pauli principle. Here we search for compact, possibly 
accurate, expressions for the eigenvalues and eigenvectors of H and to unravel 
hidden correlations among the solutions. 

We start by recasting the eigenvalue equation in the Bargmann-Fock rep- 
resentation of the hamiltonian formalism, where the odd (anti-commuting) 
Grassmann variables Xjm, replace the fermionic operators djm, a^^. For 
this purpose we introduce the even variables 

^jm = ( 1)"' Xj—m^jm y (4) 

in terms of which the operators become 

and the Hamiltonian (better, the normal kernel of) 

^ = E^^ E ^U.^rn.-Gj^n'^-^- (6) 

The index of nil-potency of the collective Grassmann variable $^ (to be 
referred to as s-quasi-boson) is fl^, i.e., 

(<l>^r = foTn>n^. (7) 
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We now search for eigenstates of n pairs of fermions in the s-quasibosons 
subspace as products of n factors, namely 

n 

M^l = U^l (8) 

k=l 

where 

L 

i/=i 

is a superposition of s-quasibosons placed in all the available levels. 

For this scope it is convenient to start from the effective Hamiltonian 

L ji, L 

n.^{^\ ^) = Y1 2e. Yl ^U.^^^^. -GJ2 ^'^^ ' (10) 

coincident with (^) in the s-quasibosons subspace spanned by the states 
Indeed while terms like A*A count the number of particles, (p*(p = A*A*AA 
counts the number of pairs. The eigenvalue equation then reads 



\ u=l mu=-iu 



L 3ii 



X 



where 



^^Pj-E E J = KV^n($*) , (11) 

^1=1 m^=-jp 



[rfA*'dA']^n n d>^*^,mJK,^. ■ (12) 



By expanding the exponentials in ([TT|), only the even powers, hence only the 
if variables, survive. Thus (pi]) can be rewritten as 



[ W]^efr(<^*, ^')M{^* + <^*', = EnM^l , (13) 
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where 



U=l m^=l/2 



(14) 



The integrals over the yp's, relevant for dealing with eq. (|TT]), are listed in 0, 
|[. One gets 



1 



fc=l 



2e^ - r]k 



(15) 



the T^fc being the solutions of the non-linear system 



n 

n, - n, n 



2e^ -Tjk j^Vi-Vk G 

l^k 



(16) 



In this paper we confine ourselves to the case of a single pair only. Then 
([T6| ) reduces to a single equation and the wave function reads 



^i($*) = 5^/3,<i>:. 



(17) 



Since the action of H on (|T^) is 



u=l /^=1 m^=l/2 

L L 



(18) 



U = l fJ,,U = l 

the Schrodinger equation becomes 

L 



(2e,-E)/?,-G^^]^/3, 



C = 



(19) 



which imphes 



L 

(2e, - E)^, -GJ2 ^t^f^t^ = Vz/ = 1, L . (20) 

11=1 

Since ^^^j^ ^^J.P^l does not depend on the index u, it follows that the coeffi- 
cients j3„ are 

C being a normalisation factor. Inserting (|2T|) into (^) we then get the 
well-known eigenvalue equation (referred to in the following as secular equa- 
tion) i 



u=l 



which yields L eigenvalues E^^^ (1 < fx < L), the corresponding components 
(31''^ of the wave function being given by eq. (pl|). 
Actually the normalisation of the state 

L L 

{ri.i^,) = (<i>.,$:o = y,^,\pa' = i (23) 

i/,u'=l u=l 

suggests to introduce the coefficients 

^l^^ = v^/3('^) . (24) 



It is straightforward to solve equation (|22D numerically: the solutions 
can be graphically displayed as the intersections of the Ihs of (^) with the 
straight line E = 1/G. In Fig. |l| this is done for 5 levels of a 3-dimensional 
harmonic oscillator and for two typical values of G = G/ItjjJq. In the figure 
two classes of states (labelled as /c = 0, . . . , A/" — 1) appear: the ffist one 
embodying the = state, which lies below the lowest single particle level 
for an attractive interaction and corresponds to a collective state; the other 
embodies the so-called "trapped" solutions (1 < A; < 4), which lie in between 
the single particle levels. 
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Figure 1: Graphical solution of the eq. (^) for the case of a harmonic os- 
cillator well with 5 levels. The values of G are 1/2 (dot-dashed) and 1/15 
(dotted). All quantities are in dimensionless units {E = E/huo). 



The collective and single particle character of the eigenstates is apparent 
in the histograms of Fig. ^, where the coefficients are displayed. The 
states on the left panel correspond to the coupling constant G = 1/2, those 
on the right panel to G = 1/15. From the top to the bottom the coefficients 
of the collective {k = 0) and of the trapped (1 < < 4) states are shown. 
The bars, from left to right, yield the components /i = 0, ■ ■ - 4 of the wave 
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-1 L -1 L 



Figure 2: The wave functions components, (3^ , of the collective and trapped 
states (from top to bottom) for G = 1/2 (left panel) and G = 1/15 (right 
panel) in the same case as in Fig. |T[ 
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functions 



M-l ~n{k) 




(25) 



In the G = 1/2 case all the l]"^^^ are sizable, reflecting the high degree of 
collectivity of the state, although the weight of the s-quasibosons living in 
single particle states with higher degeneracy, and hence more distant from 
the collective state, tends to be greater. For G = 1/15 the force is too weak 
to organise any collective motion. 

Concerning the eigenstates of the trapped solutions when G is large, the 
k=l and 2 states approach the lower unperturbed levels, and so do their wave 
functions: in fact in the left panel one sees that the dominant component of 
the latter corresponds to h = k — 1. Instead the energies of the k=3 and 
4 states occur close to the middle of two unperturbed energies: hence their 
wave functions display two dominant components. For small G the energies 
of the trapped solutions remain close to the unperturbed energies, hence their 
wave function almost coincide with the unperturbed h = k state. 

3 The collective solution 

We first consider the collective solution in the extreme situation where all 
the levels coalesce. Then eq. ( p2D reduces to 



(with Q = and = e Wu = 1, . . . , L), entailing 



Thus in this limit only the collective state survives. Clearly (|27|) remains 
a good approximation when the degeneracies of the L levels are lifted only 
if their spread in energy is small with respect to QG. An accurate analytic 
solution when the distribution of the levels is arbitrary will be derived in the 
next subsection in both the large and small G limits. 



E = 2e-nG . 



(27) 
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3.1 The 'strong coupling' limit 



For large G it is convenient to recast eq. (|22|) as 

0. I ^ n„ 1 



EQjjj ^ 1 \ ^ 
Of — P~ R — Op 



2e^ - E E-2e^ . o - e G ' (28) 

1/ 1/ i — z— 

E - 2e 

and expand in the parameter 2{ejj — e)/{E — 2e). Defining 

e = ^^^^ (29) 

in leading order we get 

Eo = 2e-nG , (30) 



which coincides with the degenerate case value (pT]). Importantly, owing to 
the definition (^), the next-to-leading order correction vanishes. 
To proceed further we rewrite (^) as 

n=0 ^ ' 71=0 ^ ' 

where the generalised moments of the distribution of single particle levels e,y, 
namely 

M(") = a"m(") = ^Y1 - ' (32) 

V 

have been introduced together with the expansion parameter 
and the variance 



^Y.^u{e.-^y (34) 



of the level distribution. The 'strong coupling' regime then corresponds to 
a -C 1. Note that the second moment is the square of the variance a, whereas 
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the third, M^^^ = a^'j, and fourth, M^'^^ = (c + 3)cr^, moments are related to 
the skewness 7 and to the kurtosis c of the distribution, respectively. 

In a perturbative scheme, setting E^^^ = E^^~^^ + 5 and linearising in 5, 
we get for the collective energy 

= -l^a^ + ^a'-{l + c)a' + 0{a') , (35) 

or, equivalently, 

an expression vahd for a <^ 1 and independent of the single particle energies 
distribution. 



3.2 The 'weak coupling' limit 

When a > 1 the collectivity is very weak (see Fig. |^, right panel) and in (^8|) 
the term u = 1 dominates the sum. Thus by separating the latter we rewrite 
([28| ) in the form 

2(e^ - ei) -{E- 2ei) ^ 2ei - E ^G' ' 



v=2 



Now, for G ^ 0, the quantity 2ex — E is of the order of G: hence the last term 
in the Ihs of (^) dominates. Thus the leading contribution to the energy is 

Eq = 2ei - Gfii . (38) 

Factorising then 2(e,y — ci) in the denominator of (^), expanding in powers 
of {E — 2ei)/2{eu — ci) and introducing the inverse generalised moments 

Af<-) = = _L-V^^ (39) 
a' fi- fit ^ (c„ - e, ' ' 

we get 



E - 2ei 



(40) 



Proceeding then as in section 3A, we next expand in powers of /5 = m*- ^^a', 
being 1/a' = G{il — i7i)/2(7. At the order 0^ we obtain 



E - 2ei 



-<(l + /5 + r 1- 



(-2) 



1^1 



1-3 



1^1 m(-3) 



(41) 



+ 



(m(-i))2 - (m(-i))3 V^-^i 



1^1 



where a' does not coincide with 1/a because the lowest unperturbed state has 
been separated out. Note that not only where the 'strong coupling' expansion 
fails the 'weak coupling' one holds vahd, but also an overlap region appears 
to exist where both expansions yield quite accurate results, being at the same 
time a < 1 and /5 < 1. 



3.3 The Euler-McLaurin approximation 

An interesting approach to the pairing problem is offered by the Euler- 
McLaurin formula It enables us to replace the sum with an integral, 



v=a 



The formula (|4^ ) holds valid if f{y\ v to be viewed as a complex variable, is 
analytic in the strip a < Keu < b {a and b being integers). Neglecting the last 
term in the rhs, we shall illustrate the use of ( ^2]) in the case of the harmonic 
oscillator well, whose levels we label with an index = 0, 1, ■ ■ ■ A/" — 1. For 
this potential, from the well-known unperturbed energies 

Cfc = ekfu^o = {k + 3/2)huJo (43) 

and associated pair degeneracies 

fifc = (fc + l)(fc + 2)/2 , (44) 

^The function R was first fixed by Euler to be /^(u — [u] — 1/2) f'{u)du. McLaurin 
provided more accurate expressions for it. 
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the total degeneracy Q, the average energy e and the variance a are found 
to be 



f]= iAr(A/'+l)(Ar + 2) 
6 



(45) 



and 



respectively. The above entail 



(46) 



(47) 



(48) 



Incidentally, from (|48|) the coincidence of ( ^6l) at the order with the finding 
of ref. PI follows. 

Before exploiting (^) we recast (|22D using the digamma functions and 
the dimensionless variables G = G/huo and E = E/huo] we get 



a 



(h.o.) 



Shuo /3 v/(A/'- l)(A/' + 3) 



G V 5AfiAf+l){Af + 2) ■ 



N=0 



iN+l){N + 2) 
2N + 3-E 



U{2 + E + Af) E^-1 



Af{2 + E + J\f) 



3-E 



3-E 



+ 



E^-1 



l + 2Af-E 



2 

G ■ 
(49) 



In the 'strong coupling' limit the collective solution is strongly pushed 
down in energy, hence the use of the asymptotic formula 



^ oo 



2k 



2z ^ 2kz'^^ 

k=l 



(50) 
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the being the Bernoulh numbers, is appropriate. It is then remarkable 



that by inserting the ip as given by the first two terms on the rhs of (|50D 
into ( ^9[ ) one obtains the same expression provided by the Euler-McLaurin 
formula for the harmonic oscillator, namely 



E-McL 



U{2 + E + U) 



+ 



1 



E 27V + 1 - E 



+ In 



2M+1-E 
3-E 



(51) 



In the tables |l] and ^ we refer to the eigenvalues obtained using (|5l|) as ML. 
Furthermore, being the term z~'^ absent in ([50|), to keep in the expansion (50) 
also the k = 1 terms yields excellent results. We refer to this approximation 
as ML2. 



3.4 Numerical results for the collective energy 

In this subsection we present our predictions for the collective energy, essen- 
tially in the harmonic oscillator case, for (5=0.1 and 0.2 (a realistic estimate 
for atomic nuclei with mass number A > 100). Our results, corresponding to 
eqs. (|35D , (0) and to the solution eq. ( ^91) with given by (0) , are shown 
in tables [l|, |^ and compared with the exact solutions for some values of A/". 



AT 


-^cxact 


0(«2) 




a 






/5 




EuL2 


2 


2.883 






2.17 


2.885 


2.883 


0.15 


2.927 


2.822 


3 


2.860 






1.34 


2.87 


2.860 


0.3 


2.909 


2.800 


4 


2.820 


3.925 


1.479 


0.89 


2.853 


2.825 


0.47 


2.877 


2.763 


6 


2.529 


3.695 


2.930 


0.46 


2.814 


2.681 


0.86 


2.590 


2.500 


8 


0.1361 


0.5375 


0.2161 


0.28 






1.35 


0.1338 


0.1362 


10 


-6.477 


-6.298 


-6.461 


0.19 






1.95 


-6.459 


-6.477 


12 


-17.68 


-17.58 


-17.68 


0.14 






2.63 


-17.64 


-17.68 



Table 1: Comparison between the exact and the approximate collective en- 
ergy for the harmonic oscillator case, for various values of A/" and G = 0.1. 
The energies are in units of hujQ. Columns 3 and 4: strong coupling; columns 
6 and 7: weak coupling; columns 9 and 10: Euler-McLaurin approximations. 
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u 




0{a') 


0{a') 


a 


0{(3) 




(3 


EuL 


EmL2 


2 


2.728 






1.1 


2.74 


2.727 


0.3 


2.796 


2.677 


3 


2.583 


3.1 


2.185 


0.67 


2.68 


2.589 


0.6 


2.653 


2.547 


4 


2.182 


2.713 


2.243 


0.44 


2.613 


2.318 


0.93 


2.224 


2.171 


6 


-1.480 


-1.303 


-1.461 


0.23 






1.73 


-1.478 


-1.480 


8 


-11.06 


-10.98 


-11.05 


0.14 






2.71 


-11.01 


-11.06 


10 


-27.94 


-27.90 


-27.94 


0.10 






3.88 


-27.85 


-27.94 


12 


-53.66 


-53.64 


-53.66 


0.07 






5.25 


-53.53 


-53.66 



Table 2: Same as table [T| for G = 0.2. 




0.1 0.2 0.3 

G 



Figure 3: The collective energy of ^^°P6 as a function of the pairing strength 
G (in MeV). Solid hne: exact numerical solution; dotted: degenerate solution 
(order dot-dashed: the approximation (|35|) at the order a^. 
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The approximate energies ( P^D are in good agreement with the exact ones 
for G > 0.2. Instead for smaller G and 4 < A/" < 8 the collective energy in 
the 'weak coupling' limit, namely eq. (|4TD, is in better touch with the exact 
values. The existence of a region where both the 'strong' and the 'weak' 
approximations hold valid, is apparent. We also quote in the tables |l|, |^ the 
results obtained with the McLaurin approximation in first and second order, 
as previously discussed. 

To address a situation of physical interest, we consider the case of 2iop6, 
where one pair of neutrons lives in the = 6 shell, set up with the s.p. 
levels 4/2, 5-7/2, Si/2, 4/2, ji5/2, Hi/2 and 5-9/2 of energies -1.5, -1.6, -2.0, 
-2.3, -2.4, -3.1 and -3.9 MeV, respectively. The total degeneracy fl, the 
average energy value e, see (^^, and the variance a then turn out to be 29, 
-2.61 and 0.77 MeV, respectively. Since for lead G ~ 0.1 MeV we obtain 
a = 0.53 and (3 = 0.89. Hence we are in the 'strong coupling' limit. 

In Fig. 1^ we compare the collective energy (|35|) for ^^°P6, at order a° 
and at order a^, with the exact numerical solution as functions of G: we see 
that, in the physical range of G, an excellent accord with the exact result is 
obtained when the term is included. 



4 Trapped solutions 

4.1 The general case 

The evaluation of the trapped eigenvalues can be easily performed numeri- 
cally. This approach however hides the interesting pattern the trapped so- 
lutions display, already suggested in where they were associated with a 
quantum number and a parabolic behaviour in the latter was shown to occur 
for a harmonic oscillator well. 

To shed light on this aspect of the pairing problem we first observe that, 
irrespectively of the strength of the interaction, 2ej,_i < E^'^^ < 2e,^, E^^^ 
being the energies of the trapped solutions {y = 2, . . . L). 

Next we define a new variable z'^"'^ according to 

E(") = 2e,_i + 22('^)(e, -e,_i) (52) 
(clearly z'^^'> G (0, 1)) and select from the secular equation the terms associ- 
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ated with the poles in 2e^_i and 26^,, thus getting 
with 



n, 



2(e^ - ej,_i) 



G 



\ 2{e^ — e^_i 



(53) 



(54) 



We have thus isolated in the Ihs of (|53D the poles trapping the solution 
in the interval (0, 1). Now we approximate the rhs. To this purpose we 
expand in powers of 2^^) (/i ^ z/,i/- 1): 



E 



n=o / - e,,_i 



n+l 



E 



Cj/— 1 , 



„=o / e,, - 



n+l 



if /X > z/ 



(55) 



if ^ < u — 1 . 



Of course (0), truncated at the order m, provides a polynomial approxima- 
tion for the function Lp. 

Next we let the discrete variable v become continuous by means of the 
Euler-McLaurin formula (|4^) . We thus obtain for ip the expression 



u-2 



E-McL 



(z/,2;) = -C(l,2;) + -C(z/-2,z)+ / C{^^,z)d^ 



2[e{u) - e{u - I)] 
G 



C(/i, z)dij, 



(56) 
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which is analytic in the variable z/, if and Vty have been analytically ex- 
tended. 

Approximate expressions for z'^^'' then follow from the continuous version 
of (0), namely 

^^ + ^ = ^™(^.). (57) 

which implicitly defines function of the complex variable z/. The 

simplest approximation corresponds to replace the rhs with (/9^~'^'^'"(z/, 0). 
One gets 

^E-McL(^^0) + fi(z/-l)+r](z/) 



2^E-McL(^^0) 

^[(^E-McL(^^ 0) + VL{u - 1) + ^](z/)]2 - 4r](z/ - l)y?E-McL(^^ Q) 



2(^E-McL(^^0) 

(58) 

which can be further simplified via a parabolic expansion around, e.g., the 
middle point v = L/2, namely 



z,{v) ^ z,{L/2) + (^^ - 2 J + 2 1,^ " 2 J ^'"^^''^^ ' ^^^^ 

To further proceed a specific model for VL{i>) and e(z/) should be chosen, 
which we shall do in the next subsection. Here we display the numerical 
results for the case of ^^^Pfo already discussed in subsec. |3]^. In tab. || we 
display both zoi^u) as given by (|^) and ziiy) (namely the ziy) obtained by 
truncating (|55D at first order and then by solving the corresponding equation 
by successive linearisations) together with the exact numerical solutions. The 
corresponding energies, see (|52D , are also reported. 



4.2 The harmonic oscillator case 

In this subsection the functions f2(z/) and e{v) are those of the harmonic 
oscillator. As in sec. p.3| , we label the solutions with the index 1 < k < 
J\f — I (the value k = corresponds to the collective solution). 

Using the dimensionless variables E, G and z = z/Tiujq we rewrite ( ^31) in 
the Euler-McLaurin approximation (^6]) obtaining 



fi(A;-l) , n{k) . 



~z{k) z{k) 



^^--'^\krzm (60) 
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V 


^cxact 


^0 


z\ 




Eo 


El 


2 


0.560 


0.621 


0.577 


-6.903 


-6.807 


-6.876 


3 


0.555 


0.605 


0.577 


-5.423 


-5.353 


-5.392 


4 


0.760 


0.762 


0.760 


-4.648 


-4.648 


-4.648 


5 


0.929 


0.934 


0.931 


-4.043 


-4.040 


-4.041 


6 


0.730 


0.749 


0.748 


-3.416 


-3.401 


-3.402 


7 


0.783 


0.782 


0.783 


-3.043 


-3.043 


-3.043 



Table 3: Exact and approximate energies for the trapped levels in the case 
of ^^°P6. The i?'s are in MeV, the z's are pure numbers. 



with 



E-McL 



(k,z) = -t- ^(4A; + 25 + 3) + ]{Af - l){Af + 2k + 22 + 3) 



+ 



k{l-k) {k + 2){k + 3) 
4(z + 1) 4(z-2) 

1 MjM+l) 

2{k + 2-l)~ A{k + 2~M) 

{2 + l){k + 2-^) 



(61) 



+ -{k + 2){k + 2 + l)\og 



{2-2){k + z-l) 



The Euler-McLaurin approximation is not only valid quantitatively, but 
also provides the key for studying analytically limiting cases (like TV — >■ 
oo), which are helpful in shedding hght on the properties of the trapped 
eigenvalues. Indeed consider eq. (^): on the Ihs it displays the usual two 
poles in the complex plane of z, at and 1, respectively, whereas the rhs has 
poles at z=-l and 2. Equivalently, in the variable E^''\ the poles in the range 
from ^^=-2) to ^(^=+1) are explicitly kept, while the others are simulated by 
two cuts. The graphical solution is shown in Fig. ^ where the quantity 



Q{k-1] 



n{k) 

z{k) - 1 



E-McL 



Gr 



(62) 



(as a function of E) is plotted. The eigenvalues correspond to the intersec- 
tions with the straight line 2/G. Four solutions are found: the collective one 
is of no interest here, the other three correspond to the indices k — 1, k and 
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Figure 4: Graphical solution of the secular equation in the Euler-McLaurin 
approximation. Only the part of the curves outside the cuts (displayed as a 
black area) is drawn. 



k + 1. We shall only keep the middle solution, the most accurate, as it lies 
far from both the left and the right cut. 

The results of the Euler-McLaurin approximation for the trapped energies 
are shown in Figs. ^ and 0, solid line. 

We now explore on the basis of ( |60D the A/" ^ oo limit. For this purpose 
we introduce the variable 

A = A. (63) 

Clearly, when Af is large, A G (0, 1). 

Also we express the coupling constant G in terms of a according to 
Then eq. (^) becomes 



1 



1 



+ 



1 



^ + log. 



1 2(5+1) 2(5-2) 
1-A 1- 



-log 



5 + 1 



fa 



A 



2A2 



(64) 
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the 5-dependence appearing only in the Ihs, referred to as 0(5), and the A 
dependence only in the rhs, referred to as ^oo(A). 

Eq. ( p^ ) is easily solved numerically and the results are displayed in Fig. ^ 
for different values of a. 




Figure 5: Solution of the eigenvalue equation in the limit A/" ^ oo in the 
Euler-McLaurin approximation for different values of a. Solid line: a = 0, 
dashed line: a = 0.2, dotted line: a = 0.35, dash-dotted line: a = 0.5, 
long-dashed line: a = 0.8, dot-dot-dashed line: a = 1.5. The exact results 
are not displayed as they almost coincide with the ones in the figure. 



Of significance is the a ^ case, which still carries the fingerprints 
of the harmonic oscillator and the a ^ oo one, which corresponds to the 
straight line z = 1. It is of great interest to follow the behaviour with a 
of the curves. The sum of the eigenvalues, or, better, of the 5(A), offers a 
guidance to globally follow this evolution. Actually this "sum rule" for the 
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trapped solutions can be exactly computed in the two hmiting cases a = 
and a = oo. Indeed from the Viete equations one has 

S(« = 0) = ;^E^(^) = i ^^^^ 

k=l 

and S(a = oo) = 1, respectively: these values set the limits for the area 
under the curves of Fig. ^. 

On the other hand the behaviour of the curves themselves is ruled both 
by the function ^ooW, which goes to — oo when A — > 1, whereas for A ^ 

UA) — 1+"" fo^«<«- with a„ = ^./?^ 0.29, (66) 
^ A^o |-oo fora>«„ " 8V5 ' ^ ^ 

and by the monotonic decrease of (j){z) in the interval (0,1), which varies 
within the limits 



Thus, when A 



\im(f){z) = +00 lira (f){z) = — oo . (67) 

2— >0 5^1 



5(A) - 1 + (68) 
41og(l - A) 



no longer depends upon a and slowly approaches 1. 

For A — >• two cases occur: if a < then z{X) ^ 0, if a > instead 
z{X) 1. Thus a transition occurs at indeed the exact eigenvalue 
arises by perturbing the {k — l)-th free one when a < acr (strong coupling 
regime), according to (|53D , and by perturbing the k-th one when a > acr 
(weak coupling regime) according to 

E(^) = 2e,-2z(^)(e,-e,_i) . (69) 

This is strikingly illustrated in Fig. |^ that shows that the almost parabolic 
behaviour of z{X) for small a is strongly distorted for a ~ for small A. 
For larger a a smoother behaviour is recovered. In particular in Fig. ^ a 
marked minimum is seen to develop for a above, but close to, the critical 
value. When M is finite, for A < l/M, there are no eigenvalues, since z lives 
on a discrete set of points: hence is ill-defined. 

A deeper insight of the above findings is offered by the following com- 
ments: 
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1. the pairing interaction is of finite range, therefore a pair trapped by 
highly excited harmonic oscillator states has the two partners extremely 
de-localised, hence unaffected by the interaction. This explains why all 
the curves in Fig. |^ coalesce to 1 when A — *■ 1, no matter what the value 
of a is. On the other hand all the curves (up to a critical value of a) 
converge to 5(0) = 0, refiecting the pressure exercised by the infinite 
number of the high-lying, large degeneracy, levels on the low-lying, 
low-degeneracy, ones. 

2. The eigenvalues obey a "sum rule" , whose value grows with a from 1/4 
to 1 , hence they must grow as well: since the action of the pairing force 
is gauged by the product GVLk^ at some critical value of a (hence for 
G sufficiently small) the system prefers to fulfil the sum rule by lifting 
the lowest eigenvalues (corresponding to the lowest degeneracies) to the 
unperturbed values. Indeed (see next section) when the degeneracy of 
the harmonic oscillator levels is artificially made to decrease, then no 
transition occurs. 

In concluding this section we turn to the problem of finding analytical 
expressions for the trapped eigenvalues. In this connection we first notice 
that, remarkably, the large M limit is still a good approximation to the 
exact solutions even for A/" = 5, as shown in Fig. ^ While in this limit the 
complexity of the Ihs of eq. (0) is still such to render difficult the finding of 
simple analytical approximations, yet its structure is dominated by the poles 
in and 1 that trap the solution. 

It is thus reasonable to approximate 0(5) with 

which retains the pole structure and reproduces (j){z) together with its two 
first derivatives at 5 = 1/2. Hence, from the equation (|6lD, 0^pp''(z) = .^oo(A), 
it follows 

= i + il - ^MlM , (71) 

displayed as a dotted line in Figs. |^ and |7|. 
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Figure 6: Exact and approximate values for the trapped solutions in the 
A/" = 5 case at different values of a: diamonds: a = 0.1, crosses: a = 0.4, 
stars: a = 1. Solid lines: Euler-McLaurin approximation, dashed lines: the 
same but in the limit A/" — >• cxo, dotted lines: approx. eq. (|7TD , dash- dotted 
lines: parabolic approximation. 



By expanding in A (say around 1/2) and in a up to the second order a 
parabolic expression for function of A is obtained, also shown in Figs. 

and 0. 

In these figures the Euler-McLaurin approximation in the case J\f —>■ oo is 
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seen to be remarkably stable and in good accord with the exact eigenvalues 
when Af is reduced to finite values for a < acr (but not for a > acr)- Further- 
more the predictions of the simple expression eq. (^) almost superimpose to 
the Euler-McLaurin ones. 



z(k) 



0.8 



0.6 



0.4 



0.2 




10 



J I I I I L 



15 20 
k 



Figure 7: Same as in Fig. M but with Af = 15 
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4.3 Variations on the harmonic oscillator 



In this subsection we explore whether the transition previously discussed 
occurs only in the case of the harmonic oscillator or is more general. Ac- 
cordingly, we take M large and leave unchanged the regular arrangement of 
the harmonic oscillator eigenvalues, but vary, admittedly artificially, their 
degeneracy according to the prescriptions 



(a) 
k 



(Ar-A;)^~Ar^(l-A)' 



(72) 
(73) 



with 7 > 0. From the above 



'(7 + 1)= 



(7 + 3) (7 + 2)G^^^ ■ 



(74) 



Note that for 7=2 (|7D yields a 



17=2 



Q,('!-o.)/2 in the large M limit. 



Now from eq. (|56| ) one gets for the function ip^ '^'^'^^{k, z) in the M ^ 00 
limit and in the case (a) 



E-McL 



{k,z) 



AT^oo 2 

+ Ar^<'p(A,7) + A^lo 



3 + 7 



z + 1 ' Z-2J ■ ■ V (1 + 7)3 
(l + z)(l-A) 



;2 + 7)a 



(a) 



(2 - z)X 



(75) 



(a similar expression holds for case (b)), where 



(76) 



is an analytic (hence well behaved) function of A and goes to I/7 when A — > 0. 
Accordingly the secular equation can be recast in the form 



log 



1-A 



1 

A^ 



P(A,7)- 



7 + 3 
(7 + 1)^ 



■(7 + 2)a(") 



(77) 
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in case (a) and 



log 



1-A 
A 



A)' 



P(l-A,7) + 



7 + 3 
(7 + iy 



■(7 + 2)a 



(78) 



in case (b), the Ihs exactly coinciding with the one of eq. 

In the limits A — >■ (case (a)) and A — > 1 (case (b)), the leading terms of 
(|73) and (^ are 



1 

A^ 



7 + 3 
(7 + 1)^ 



■(7 + 2)«('^) 



(79) 



and 



(1 - A)' 




7 + 3 
(7 + 1)^ 



■(7 + 2)a 



(80) 



respectively. Remarkably z{X) in the A — > limit goes to or 1, according 
to the sign of ([79|) ( ), which is set by a. The critical value is 



a 



(a) 



'(7 + 1)' 



7 + 3 7(7 + 2) 



(81) 



clearly behaving as I/7 when 7 — 0. On the contrary, in the case (b) in the 
A — >■ 1 limit z{X) always tends to 1, since (|80[) never changes sign. 

Thus a transition occurs only if the degeneracy is growing with k and the 
'strong coupling' domain becomes wider («„ increases) as 7 approaches zero: 
here the transition disappears. If the degeneracy decreases with k (case b), 
no transition exists. 

The eigenvalues corresponding to the degeneracies (a) and (b) for a = 
are displayed in Fig. ||. 



5 Conclusions 

In this paper we have once more addressed the pairing problem with the 
scopes of finding approximate, but analytic, solutions for the eigenvalues and 
eigenvectors and of disclosing a surprising pattern displayed by the former. 
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Figure 8: The functions z{X) for different degeneracies and a = 0. Solid line: 
case a), 7 = 2; dashed hne: case a), 7 = 1; dotted hne: case a) and b) 7 = 0; 
dash-dotted line: case b), 7 = 1, dash-dotted-dotted line: case b), 7 = 2. 



To pave the way to the general problem with any number of pairs and 
levels we have addressed the simplified situation with only one pair living in 
a set of Af single particle levels of a potential well. 

For the collective mode we have found, in the strong and weak coupling 
regimes, expansions relating its energy not only to the strength of the pairing 
interaction, but as well to the parameters characterising the distribution of 
the single particle levels. Thus our analysis shows that different potential 
wells could lead to the same energy for the mode, which is accurately fixed 
by only a few parameters of the levels distribution and not by the precise 
energy of each of them. In particular the variance of the levels distribution 
appears crucial for the development of a collective mode. 

Concerning the energies of the trapped solutions (when these are known. 
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so are the related eigenfunctions) our search for a simple formula for their 
description has in part been prompted by the recent finding concerning 
the integrability of the pairing Hamiltonian, both at the quantum and at 
the classical level. Indeed it is reasonable to expect the existence of rela- 
tively simple analytic formulas for the eigenvalues of a Hamiltonian if the 
corresponding classical motion is not chaotic. 

For this scope the Euler-McLaurin approximation, conveniently exploited, 
has been invaluable. We started from the extreme case of an infinite number 
of single particle levels of an harmonic oscillator and classified the trapped 
solutions in terms of a quantum number A, varying between and 1. For 
sufficiently large values of the pairing strength G the displacements of the 
trapped eigenvalues from the unperturbed solutions always start from and 
monotonically grow to 1, reached when A=l. The first of these findings 
relates to the "pressure" exercised by the high-lying eigenvalues (associated 
with very large degeneracy) on the low-lying ones, where the degeneracy is 
low. The second one instead relates to the delocalisation of the partners 
of the pair in very high-lying harmonic oscillator levels. Remember indeed 
that the pairing interaction is meant to simulate the short-range part of 
the nucleon-nucleon force: hence it is incapable of correlating two fermions 
lying far away from each other, provided no quasi-bound is generated by the 
interaction, as is the case for the trapped states. 

In this connection it is worth reminding that the classical limit is achieved 
by letting the degeneracy of the single particle levels become very large 0. 
Thus the coincidence of all the eigenvalues, for any G, in A=l also refiects 
the evolution from quantum to classical mechanics of our system. 

Another aspect of significance of our work concerns the transition of the 
behaviour of the eigenvalues from one regime to another. Indeed we proved 
that a critical strength of the pairing force exists such that for weaker in- 
teractions the behaviour of the trapped eigenvalues versus A ceases to be 
monotonic. Actually for all G < Gcru the displacements of the eigenval- 
ues from the unperturbed energies start from 1 rather than from 0. This 
occurrence might be understood on the basis of the sum rule the trapped 
solutions should fulfil and of the nature of the pairing force. In fact at some 
point the system prefers to obey the sum rule by lifting the lower eigenvalues 
(associated with low degeneracies) to values close to 1 (where the levels are 
associated with very large degeneracies). 

Finally, from our analysis it emerges that the eigenvalues obtained in 
the —>■ oo are very robust with respect to variations of A/", when G is 
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large: indeed they keep their vahdity even for values of M as small as 5. 
Moreover, for G > Gcrit, we have found that the trapped eigenvalues indeed 
lend themselves to simple analytical expressions, even to a parabolic one, as 
hinted in ref. 0. 

We are presently investigating the statistical fluctuations of the trapped 
eigenvalues, which should reflect the integrability or, equivalently, the ab- 
sence of chaotic motion associated to the pairing Hamiltonian. 
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